Pinning of slidding collective charge state in a ID attractive fermion system 
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We investigate an interacting fermion model with boundary potential by using Bethe ansatz 
method. The ground state properties of the system and the boundary effect are discussed. It 
is found that attractive boundary potential leads to the boundary bound state. An interesting 
phenomenon is that the slidding collective charges in a periodic system, which is formed due to 
the attractive interaction among the fermions, will be pinned around the boundary, as long as the 
negative boundary potential is strong enough. 
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I. INTRODUCTION 

There has been extensive interest in the investigation of low-dimensional correlated fermion systems in the recent 
years. It is well known that the perturbation techniques are not very effective tools in dealing with the one-dimensional 
systems. The ID systems, although they are somewhat artificial, can give some valuable information on the role of 
correlation effects in higher dimensions and thus are theoretically meaningful. Many approaches, such as bosoniza- 
tion and renormalization group techniques, have been successfully applied in this field and some fruitful results are 
obtained. Usually, studies on the exactly solvable models can provide an exactly theoretical understanding to the 
more complicated interacting systems. Thus the quantum integrable models, which have a long history in itself, again 
attract much attention. 

Some recent works Jl],|| show that the behavior of a single impurity in the ID quantum systems is rather different 
from that in a Fermi liquid. Kane and Fisher || first investigated a ID repulsive interacting system in the presence of 
a potential barrier and pointed out that it corresponds to a chain disconnected at the barrier site at low energy scales. 
This can be effectively described by the open boundary condition and is well investigated by the boundary conformal 
field theory Q . The open boundary problem is also studied by using the Bethe ansatz method (5|-0| . Attempting to 
understand the effect of boundary potential to an interacting fermion system, we study a spinless interacting model 
with boundary potential by the Bethe ansatz method. The results show that there exists a critical value for the 
attractive boundary potential. When it goes through the critical value, the strong attractive boundary potential leads 
to a boundary bound state. Moreover, the attractive interaction among the fermions also has a critical value. When 
the interaction is strong enough, a collective slidding charge state forms. The attractive boundary potential will 
make the whole collective state to be pinned down around the boundary barrier site. This phenomenon provides us 
a meaningful example how a local potential influences the global properties in a strongly correlated fermion system. 

II. THE MODEL AND BETHE ANSATZ 

Consider a ID lattice with L sites and N particles. The Hamiltonian reads 

H = T,j=i[-' t ( a i a j+i + a i+i a i) + UrijUj+i] +pxnx +p L n L , (1) 

where aj(a^) are fermion annihilation (creation) operators and t, U, p\, Pl are hopping amplitude, coupling and 
boundary potential constants respectively. The hopping parameter is defined to be positive (t > 0). As a matter 
of fact, this model is equivalent to the XX Z Heisenburg model ||||by a Jordan- Wigner transformation |l(J. For 
simplicity, we take pl = oo, which means that the particles arriving at the right boundary are completely reflected 
by a infinite high wall, and only consider the effect of the left boundary. To make the discussion clearer, we rederive 
some known results about this model in this section. 
Taking the eigenstate as 

I = ( Pi x ^'-- x N)at 1 •••<4 JV | 0), 

X\---Xn 
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we get the eigenequation 



-tY,v{xi-Xk ± 1, -x N ) + 8 xl+ i )Xk ip(xi ■ --x N ) + [pi Y^Hxk - 1) 

k k<l k 

+PN S ( x k - L)]ip(xi ■ ■ ■ x N ) = Eg(xi ■ ■ ■ x N ). 

k 



(2) 



Considering the two-particle state, in the region of x\ < x 2l we take the wavefunction as Bethe type 

ip{xi,x 2 ) = Eri,r2^i2 (ri , r 2 ) cxp(ir 1 k 1 x 1 + ir 2 k 2 x 2 ) + A 21 (r 2 , n) cxp(ir 2 k 2 x 1 + ir x k 2 x 2 )\ (3) 

where ri i2 = ±1 indicates that the particles move toward right or left. Far away from the boundary (2 < x\, x 2 < L), 
when the two particles are not neighboring (xi ^ x 2 — 1), it is easy to get the eigenvalue from the eigenequation 

E = -2t (cosfci + cosfc 2 ) . (4) 

When the two particles are neighboring [x\ = x 2 — 1), solving the eigenvalue equation, it is readily to obtain the 
following relation 

-42i {r 2 ,n) = S(r 1 k 1 ,r 2 k 2 )A 12 (n,r 2 ) (5) 

with 

qfnk-i r-.k-A — pW(.rik u r2k 2 ) — l+cxpi(rik 1 +r 2 k 2 ) + U/tcxpir 2 k2 (f .\ 
D{riKl,r 2 K 2 ) — e - i+c^pi^^+rlk^ + U/tcxpink! \°) 

Comparing to the periodic boundary case, the wavefunction with open boundary includes some coefficients cor- 
responding to backscattering (r = —1) waves. In addition, the boundary conditions give some limitation to these 
coefficients. When a particle is located at the left boundary (xi = 1) , the eigenvalue equation implies the relation 
A12 (-,r 2 ) = S L (ki)A 12 (+,r 2 ) or A 21 (-,n) = S L (k 2 )A 21 (+,n) with 

SL(k i ) = - e ~2+ Pl L t e ak *. (7) 

e ik t -\-p l /t 

For a particle at the right boundary (x 2 = L) , the eigenvalue equation also implies the relation Ai 2 (n,— ) = 
S R {k 2 )A 12 (n, +) or A 21 (r 2 , -) = S R {kx)A 21 (r 2 , +) , with 

S R (h) = -e^ L . (8) 

In general, we suppose the wavefunction has the following form 

n 

<p(x lr --,x N ) = (-l) I{Q) Ap(r p )l[e"^ k ^9(x Ql <---<x Q „), (9) 
p-Q-Xp *=l 

where I{Q) is the parity of the permutation of Q. With the same procedure, a series of relations can be obtained by 
using the eigenequation (2). Similar to the two particle case, the eigenvalue is given by 

N 

E = -2t^2cosk 3 . (10) 

When two particles are neighboring, we get 

A Pl .._p j+1 p j ._.p N (ri , • • • r pj +1 7 T pj 7 ' ' ' r p N ) = S(r Pj kp. , r PJ+1 kpj +1 )A pl p jPj+1 p N (r\ , • • • r Pj , r Pj+1 , • • • r PN ) . (11) 
with S(k Pj , k Pj+1 ) given by (6). When the particle is at the boundaries, the eigen equation gives 

A> (-,•••) = S L (k Pl )A p (+,•••), A p (---,-) = S R {k PN )A p (■■-,+), (12) 
with SL(ki), S R (ki) given by (7) and (8) . From the relations (10-12) we conclude 

A Pl ... PN {kp 1 )S{k P2 , —kp x ) • • • S(kp N , kp 1 )S R (kp 1 )S(kp 1 , k PN ) ■ ■ ■ S(kp 1 , k P2 )Ap 1 p 2 ...p N . 
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For simplicity, we use A pl ... Pj ... PN to represent A P1 Pj PN (+,■■■,+,■■■ +) and A pi ... p -... PN to represent 
A P1 Pj PN (+, •••,—,•••+) in the following text. By taking p\ = j, we get 

N 

SLi-k^Snikj) e^'-^e* 9 ^*^ = 1. (13) 

Representing S^ki) = - e i6L ^e i2ki and Sn(ki) = - e l2k * L , we have 

n 

-ikj2(L-l) = TT e i6 L {-kj) e ie{kj,hi) ^(.h-kj) _ ( 14 ) 



e 



For convenience, we take t = 1 and put A = ^ in the following discussion. A > or A < corresponds to the 
repulsive or attractive interaction respectively. As is well known, A = ±1 are two critical points. | A |> 1 indicates 
the strong coupling regions while | A |< 1 indicates weak coupling region. In the next section we discuss how the 
boundary potential influences the properties of the ground state of the whole system in different parameter regions. 

III. BOUNDARY BOUND STATES 

We discuss first the gapless critical region. For | A |< 1, it's convenient to represent A = cos 77 (U = 2 cos 77) with 
77 £ (0, 7r). The repulsive or attractive interaction is decided by 77 £ (0,7r/2) or (ir/2,ir) respectively. We parameterize 
kj as kj = 4>(Xj, tt) with 4>{a,b) defined by </>(a, 6) — 2arctan(tanhacot&). It's readily to get 6(kj,ki) = 4>{\i — Xj,r)). 
Put 

= e ir+«„ _ 1 ^ sm l(r + 77) 
Pl e ir -e'"? sin|(r-J7)' 

It follows that e l9 i( fc j)+ ,fc j — — sinh(Aj — ^-)/ sinh(Aj + The Bethe ansatz equation is thus reduced to 
( svah(\j + ±ri)\ sinh(Aj - f ) _ -pj- sinh(A 3 - r\ + irj) 



sinhiXj - fry) J sinh(Aj + f ) ^ ^1 s inh(A J - rA» - ijj) ' 



and the eigenenergy is 



N ■ 2 

sm 77 



£ = 2tV[cos, ? ] (17) 

^ L ' cosh 2 A, - cos 77 v ; 

j=l 3 I 

It should be notified that the Bethe ansatz equation with open boundary condition is reflecting invariant, which means 
Xj and — Xj corresponds to the same state, so we only need to choose one in dealing with concrete problem. 

From equation (16) , we observe that Xj — (or — i~\ is a possible solution of the Bethe ansatz equation in 

the condition of | ^° f I > 1 (I Pi I > l)when L — ► 00. It can also be seen from (17) , when cosT > cos 77 , the 

energy contributed by the boundary mode is lower than that of a real mode. For example 77 6 (0,7r/2), taking 
< |(r + 77) < |, and -| < |(r - 77) < 0, the condition | pi \ > 1 is satisfied. By the limitation < T < ~ , the 
condition cosT > cos 77 is satisfied if T < 77. This implies p\ < — 1. 

So far, we have learned that repulsive or the small attractive (—1 < p\ < 0) potential does not produce the boundary 
bound state. The BA equation has no imaginary solution in the ground state. The effect of boundary potential of 
this case has been extensively investigated by many authors ||. The strong attractive (pi < — 1) potential produces 
a boundary bound state and the corresponding energy contributed by the boundary bound state is 

_ . sin 2 77 

e/j = 2t[cos77 



cos r — cos 77 ' 



with 0<r<77<-|. In this case, the solutions of BA equation to the ground state consist of a boundary imaginary 
mode and N — 1 real mode. This problem has also been discussed by many authors pd|-p^] for XX Z and Hubbard 
model. 
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For A = — 1, the interaction term is attractive. In this case, the model corresponds to the well known XXX 
ferromagnetic spin chain [l4| with fixed magnetization. Parameterizing kj = 2arccot2A :; , we have 9{kj 7 ki) — 
2arctan(Aj — A;). Putting 



fi = I3^. ( 18 ) 



the Bethe ansatz equation is given by 



The eigenvalue is given by 



N 

With periodic boundary condition, all possible string solutions of the Bethe ansatz equation are A™ k — A„ — |(n + 
1 — 2k) with fc = 1, • • • , n. Notice x^+i/i = ~ Thus a n-string carries the energy 

^ 4 _ n 

£[ 4A^, fe 2 + 1 = A S 2 +n 2 /4' 

Representing M„ as the number of n-strings, the relation N = X^^Li "nM n holds. It can be proved 

AT N ~ 1 

A 2 + 7V 2 /4 < ^ S 2 + n 2 /4' (21) 

Thus the iV string state is favorable since it contributes much lower energy than the other states. The corresponding 
ground state wavefunction describes a bound state which decays exponentially as a function of coordinate differences 
Xj — Xi | with j ^ i. 

It can be seen from (19) , A_, = |T is a solution of Bethe equations in the limit of L — > oo if \p±\ > 1. Owing to 
Xj = |T is a solution of Bethe equation, the following type of strings 

A£ = |(r + 2fc-2) (22) 

with k = 1, • • • , n are also possible solutions. The energy of an imaginary mode is 

4 4 



4(A£) 2 + 1 l-(r + 2fc-2) 2 ' 

Now, it's easy to see that when T > 1 the energy contributed by the boundary string is lower than that of a real 
mode. 

As discussed, the attractive boundary potential produces a boundary bound state as soon as it is stronger than the 
hopping amplitude. In the case of attractive interaction, the ground state of systems will be a boundary string state 
with length N in the form (23) . In virtue of the relation 

2 y I ^ L_ 

i-(r + 2fc-2) 2 r + 2n-i r-i' 



the ground state energy is given by 

N 

E = t*p( , * 2)=2i( - —)~2Nt. (23) 

^4Af) 2 + i ; v r + 2iv-i r-r v ' 



4 _ i i 



k=l 
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So far, we have shown all the particles are bounded by the attractive boundary potential in the ground state. The 
wavefunction decays exponentially away from the left boundary. In order to show it explicitly, we give an example of 
two particle case in the appendix. Generally for a n particle system, the ground state is a boundary n-strings (22), 
the corresponding wavefunction is given by 

<p{x u -x n ) = cf (A?)-* 1 / (A£p 2 • • • / (X n n r X3 , (24) 

with /(A,) = e lk] = . Therefore, with the presence of the attractive impurity potential, the slidding collective 

charges will be pined and loses its mobility. By the expression of energy (24) , we learn that the stronger the attractive 
potential is, the more quickly the wavefunction decays away from the boundary. 

For the repulsive interaction U — 2 (A = 1) , we re-parameterize kj = 2 arctan2Aj, thus 9(kj,ki) = 2 arctan(A, — Xj). 
Putting pi = f^zj, the Bethe ansatz equation is also given by (19) , but the eigenvalue is given by 

Similar to the case A = — 1, the Bethe ansatz equation has the imaginary solution Xj = (or — i^) if | p\ > 1. 
Thus AjJ = |(r + 2k — 2) is also the possible solution of BA equation. From (25) , we see if T £ (0, 1) the energy of the 
imaginary mode A 1 — i\ is lower than those of the real modes. However, the other imaginary solutions for n > 1 arc 
not favorable in energy and thus corresponds to the highly excited state. In this case, the ground state is composed 
of an imaginary mode A = with r £ (0,1) and N — 1 real mode. 

In the following, we discuss the ground state properties of strong interaction cases, | A |> 1. First, we consider the 
attractive interaction case. Owing to A < —1, we can parameterize it as A = — cosh 77 with r\ > 0. By putting 

kj = 2 arccot(tan A., coth ^77), (26) 



6(kj,ki) — 2 arctan[tan(Aj — A^) coth 77], (27) 

_ e r +" - 1 _ sinhl(r + ?? ) 
Pl ~ e r -e" " sinh±(r-77)' [2 *> 



the Bethe ansatz equation can be reduced to 



sin(Aj + ^iy) Y" * sin(Aj- - i±) = yr jL sin(A 3 - r\ + irj) 
^m{Xj-\ir,)J sin(A, +^) ^ fj. sm{Xj - rXi - irj) " 

The eigenvalue is given by 

E = 2tY[ -— V - -cosh??] (30) 

^ L cosh?7-cos2A 7 - n v ; 

j=i j 

In the attractive interaction case (A = — cosh 77), the energy contributed by a string is lower than those of the real 
modes. In the period boundary condition, all the possible types of string solutions are 

Xl k = X n a - l - V (n+l-2k) (31) 

with k — 1, • • • , n. The energy of an n-string is given by 

e„(A)=sinh?7 sinhnry ^„ ncosn7 ^ (32) 

cosh n?7 — cos 2 A 

It can be proved that 
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sinh Nn sinh nil 

< E M «; 



cosh Nn — cos 2 A cosh 7777 — cos 2 A 

71=1 

Under the periodic boundary condition, the ground state in strong attractive case is given by an N-string |p^| , |l6| . 
Similar to the case A = — 1, the ground state wavefunction decays exponentially as a function of coordinate differences 
Xj — Xi | for all j ^ i. 

The Bethe ansatz equation (29) has an imaginary solution A^ = for T > in the limit L — ► oo. This solution 
corresponds to the boundary bound state. The imaginary mode contributes the energy 

eb(«x) = smn rj _ cosh77. (33) 

y 2 ! coshr7-coshr ' v ; 

For F > 77, the imaginary mode contributes much lower energy than that of a real mode. Similar to the discussion 
for A = — 1 case, if is a solution of BA equation, the following string 

A£ = 1^ + 2(^-1)7?] (34) 
with k = 1, • • • , n, are also solutions. The energy of an imaginary mode is written as 

sinh 2 77 

cosh 77. 



fc coshr? - cosh[r + 2(k - l)rj] 

Now, it's easy to see that the energy of this state is lower than thoes of any other states. In this case, the ground 
state solution is not a slidding ./V-string as (31) but a boundary iV-string given by (34). Put r = T — (N — 1)77. In 



virtue of 

2 



Esmh 77 . sinh nr\ 
— - = smh?7 7 , 
. cosh 77 — cosh [r + (2k — n — \)r)\ cosh 7777 — cosh T 

we obtain the ground state energy as 

N 

E M 2t sinh 77 sinh iV77 

i h coshiV?7-cosh[r + (A^- 1)77] ' 

Similar to the A = — 1 case, the ground state wavefunction is given by 

V >{x l ,--- x n ) = cf (AIT* 1 / WT X2 ■ ■ ■ f (Kr X3 , (35) 

with 



f(x n )= sinh l[r+(2fc- 1)77] 
1 [ k> sinh i[r+(2fc- 3)77]' 

For A = cosh 77 (77 > 0), we introduce the notations 



kj = 2 arctan(tan Xj coth -77), (36) 

9(kj,ki) = 2 arctan[tan(Ai — Xj) coth rj\, (37) 

e r +" - 1 sinh^r + 77) 

Pi t = — = . Y7= r. (38) 

e 1 - e 1 * sinh 5 (T — 77) 

The eigenvalue is expressed as 

N sinh 2 77 

E = 2t Y [cosh 77 S V - ] (39) 

^ L ' cosh?7 - cos2A, J v 1 

j=i j 

In this case, the string solutions correspond to highly excited states. Now, we can see Xj ■ = is a possible solution 
of BA equation (29) if the condition | p\ |> 1 is satisfied. It is easy to see that the energy of the imaginary mode 
Xj = ij is smaller than those of real modes when < T < 77 (pi < —1) for the repulsive interaction (A > 1) . 
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IV. CONCLUSION 



In summary, we study an interacting spinless fermion system with boundary potential. By analyzing the Bethe 
ansatz equation and the cigcncncrgy, we discuss the boundary effect of the ground state. Our results show that, 
when the attractive potential increases to the value of the hopping amplitude, there will be a boundary bound state 
which leads to a fermion to be bounded on the first site. When the attractive potential is weaker comparing to the 
hopping amplitude, the attractive potential is not enough to resist the hopping charges. For the attractive interaction 
U < 0, there exists a critical value U = —2t (A = —1). When the attractive interaction exceeds the critical value, 
it makes all particles correlate with each other in a form similar to the "Cooper's pair". The corresponding ground 
state wavefunction decays exponentially as a function of coordinate differences | Xj — Xi | for all j ^ i. In this case, the 
boundary bound state will influence the properties of the whole system drastically. Owing to the strong correlation 
effect among the fermions, the collective charge composite is pinned down in the neighborhood of the boundary barrier 
site when the attractive potential is strong. This is much different from the common Fermi liquid picture, in which a 
single potential barrier can not change the global properties of a system. This exactly solved model, although is very 
simple, reveals a meaningful picture of a strongly correlated fermion system. 

Appendix 

For the two particles system, the wavefunction is given by (3). We note A 2 i (r 2 , r\) — S(riki,r 2 k 2 )Ai 2 (ri, r 2 ) with 

ai i i \ riAl ~ r2A2 ~ i 
S{r 1 ki,r 2 k2) 



nAi — r 2 A 2 + i 
For convenience, we represent 



/(A,) = e«> = 



il , ^3 + 2 



\ ■ 1 

A 3 ~ 2 

thus the wavefunction can be expressed as 

V (x u x 2 ) = ]T A 12 (n,r 2 ) [/ (XiY lXl f (X 2 p X2 + S(nh,r 2 k 2 )f (\ 2 p Xl f (Xip* 2 }. (40) 



Using the relation (12) and (13) , we have A j2 = S\A 12 , = S%Ai 2 , A^ = Sl{Xi)Sb(X 2 )Ai 2 . We note 

S L (Xi) 



(A 1 + i)(A 1 + ir) 



(A! - i)(Ai - |r) 

Using the Bethe ansatz equation (14) , we have. 

S R (k 2 ) = -e t2fc2i = S L (A 2 )Si2%. 

Thus the wavefunction is written as 



<P (xi,x 2 ) = A 12 [f (Ai) 1 / (X 2 p + S 12 f (A 2 ) / (Ai) ] + A T2 [f (Ar)^ 1 / (A 2 ) + S- 12 f (A 2 ) / (A 1 )^ 2 ] + 
A- 2 [f (X,r I (X 2 y X2 + S 12 .f (A2)-* 1 / (X^ 2 } + A- 2 [f (Ai)~ Xl / (X 2 y X2 + %/ (A2)-* 1 / (Xi)~ X2 ], 



(41) 



with A- 2 = S L (Xi)A 12 , A 12 = Sl{X 2 )Si 2 S 21 A 12 , A^ = S L (Xi)S L {X 2 )Si 2 S 21 A 12 . If we take Ai = |r, A 2 = i(T + 2), 
then 5i2 = 00, Sl(Xi) = 00. That means only one term does not disappear for A 2 — Ai = i. The wavefunction has 
the form 

<p(x u x 2 ) = cf (x,)-^ f {x 2 y X2 

If we take Ai = — A 2 = — |(r + 2), S\ 2 = 0, Sl{X\) = 0, the wave function is written as 

<p(x 1 ,x 2 ) = c f(x 1 y i f(x 2 f\ 

The wavefunction can also be expressed as 

<p (x u x 2 ) = An ( r 2,ri) [f {X 2 p Xl f {X^ + S-'inh, r 2 k 2 ) f (X,)^ 1 f (X 2 p X2 }. 

r\,r1 
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It can be written as the similar form of (41) with A 21 = Sl{M)-A21, A 2 j = Sj 2 Sl(Xi)S 12 1 A21, A^ = 
5l(Ai)S , l(A 2 )S'2 1 ^42i- If we take A 2 = §I\ Ai = |(r + 2), S21 = 00, Sl(A 2 ) = 00, it is reduced to 

<p( a;i ,a ;2 ) = C /(A2)- Xl /(A 1 )- X2 . 
If wc take A 2 = — |r, Ai = — |(r + 2), S21 = 0, 5l(A 2 ) = 0,the wave function is written as 

<p(x 1 ,X2) = cf{\2) Xl f(\l) X '. 

So far, we can see, the wavefunction can be always represented as the form 

<p(x u x 2 ) =cf(xi)- xi f(xiy X2 





4) 


- Xl /r- 


3. -x 2 






U- 


-1/ 
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